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Abstract 

In these lectures notes, we will discuss a particular higher order gravity theory, Love- 
lock theory, that generalises in higher dimensions than 4, general relativity. After briefly 
motivating modifications of gravity, we will introduce the theory in question and we will 
argue that it is a unique, mathematically sensible, and physically interesting extension of 
general relativity. We will see, by using the formalism of differential forms, the relation of 
Lovelock gravity to differential geometry and topology of even dimensional manifolds. We 
will then discuss a generic staticity theorem, quite similar to Birkhoff 's theorem in general 
relativity, which will give us the charged static black hole solutions. We will examine their 
asymptotic behavior, analyse their horizon structure and briefly their thermodynamics. 
For the thermodynamics we will give a geometric justification of why the usual entropy- 
area relation is broken. We will then examine the distributional matching conditions for 
Lovelock theory. We will see how induced 4 dimensional Einstein-Hilbert terms result on 
the brane geometry from the higher order Lovelock terms. With the junction conditions 
at hand, we will go back to the black hole solutions and give applications for braneworlds: 
perturbations of codimension 1 braneworlds and the exact solution for braneworld cos- 
mology as well as the determination of maximally symmetric codimension 2 braneworlds. 
In both cases, the staticity theorem evoked beforehand will give us the general solution 
for braneworld cosmology in codimension 1 and maximal symmetry warped branes of 
codimension 2. We will then end with a discussion of the simplest Kaluza-Klein reduction 
of Lovelock theory to a 4 dimensional vector-scalar-tensor theory which has the unique 
property of retaining second order field equations. We will comment briefly, the non-linear 
generalisation of Maxwell's theory and scalar-tensor theory. We will conclude by listing 
some open problems and common difficulties. 
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1 An introduction to Lovelock gravity 

1.1 Modifying general relativity 

A convenient starting point for treating modifications of gravity are the fundamental build- 
ing blocks of general relativity (GR) itself. According to Einstein's theory, gravitational 
interactions are described on a spacetime manifold by a symmetric metric tensor g endowed 
with a metric and torsion-free connection (by definition a Levi-Civita connection) that obeys 
Einstein's field equations. In component language these equations read, 

Gab + ^bare9ab = SnGTab (1) 

where the Einstein tensor Gab = Rab ~ \9abR is given with respect to the Ricci curvature 
tensor Rab and we have included A, the cosmological constant, and T^i^ the energy-momentum 
tensor. The field equations are acquired from the Einstein-Hilbert action, 

where the Langrangian is the functional 

C{M,g,V) = -2Kbare + R, (3) 
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by variation with respect to the metric g and adequate boundary conditions (see for example 
the appendix in [1]). The bare cosmological constant A is a free parameter of the theory. 

We expect Einstein's theory to break down at very high energies close to the Planck scale, 
rn?pi = Yg^, where higher order curvature terms can no longer be neglected. Theories such 
as string theory or quantum loop gravity or again models, of extra dimensions, consider or 
model the effect of such modifications. GR on the other hand is very well tested at the so- 
lar system and by binary pulsar data in the regime of weak and strong gravity respectively 
[2]. However, recent cosmological experiments, or astrophysical data, such as galactic rota- 
tion curves, or even the Pioneer anomaly, appearing just beyond solar system scales, could 
question the validity of GR even at classical scales at large enough distances. In particular, 
recent cosmological evidence, coming essentially from type la supernovae explosions [3], point 
towards an actually accelerating universe. Looking at (1) there are 3 theoretical directions 
one could pursue in order to interpret this result. Firstly, we can postulate the existence of 
an extremenly small positive cosmological constant of value, A 

now ^ 

(lO-^eF)^, fixed by the 

actual Hubble horizon size, driving the acceleration in (1). To get an idea of how tiny this 
constant is note that this minute energy scale is most closely associated to the mass scale of 
neutrinos, lQ~^eV. Hence, although such a possibility^ is the most economic of all, since we 
can fit actual multiple data with the use of a single parameter, it actually demands an enor- 
mous amount of fine-tuning. Indeed, from particle physics, the vacuum energy contributions 
to the total value of the cosmological constant are of the order of the ultraviolet cut-off we 
impose to the QFT in question. It can therefore range as far up and close to the Planck scale 
(for discussions on the cosmological constant problem and ways to explain it see [4]). The 
"big" cosmological constant problem is precicely how all these vacuum energies associated to 
the GUT, SUSY, the standard model etc are fined-tuned each time to by an exactly oppo- 
site in value bare cosmological constant A^a^e appearing in (2). The unexplained small value 
of the cosmological constant Know is then an additional two problems to add to the usual 
" big" cosmological constant problem, namely, why the cosmological constant is not cancelled 
exactly to zero and why do we observe it now. 

A second alternative explanation one can consider, is that the accelereted expansion is due 
to a cosmological fluid of as yet unknown matter, dubbed dark energy, such as a quintessence 
(scalar) field with some potential appearing in the right hand side of (2). One of the basic 
strengths of this approach is its simplicity and in some cases an interesting approach to the 
cosmological coincidence problem. Among its basic weaknesses, apart from the usual generic 
fine-tuning and stability problems to radiative corrections, is that if we sum- up the as yet 
undiscovered matter sectors of the Universe ie, dark matter and dark energy, we conclude that 
only a mere 4% percent of the actual matter that constitutes our universe in its actual state has 
been discovered in ground based accelerators! Although there exist theoretically motivated 
dark matter candidates, such as neutralinos or axions, stemming from well motivated particle 
theories, our understanding of dark energy, is rather poor. A third, far more ambitious 
alternative, that is less well studied, far more constrained, and admittedly less succesful up 
to now, is to modify the dynamics of geometry on the left hand side of (2) not only in the 
UV but also in the IR sector. This then would mean that Einstein's theory is also modified 
at large distances at the scale of the inverse Hubble scale of today as measured in a LFRW 
universe {Hq/c = 7.566 x 10~'^'^m~^). This distance scale is enormous, to get an idea if we 

^It is not an explanation until we find a precise mechanism of why it is there at all and why now. 
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consider as our unit the distance of the earth to the sun (1 AU) we get ^ a present horizon 
distance of W^^AUl 

Next question is how do we modify gravity consistently? One can consider three basic 
types of modification which at the end of the day are not completely unrelated. Indeed, 
we can include additional fields or degrees of freedom, for example, scalar or vector, see for 
example [5]), [6]), we can enlarge the parameter space where the theory evolves, for example 
the number of dimensions, the geometric connection in question (we include torsion etc) or 
again we can generalise the field equations. 

In all cases, it is very important to fix basic consistency requirements for the modified 
gravity theory. To fix the discussion we can ask for three basic requirements: first we would 
like that the theory under consideration be consistent theoretically, for example we ask for 
sensible vacua of maximal symmetry, such as Minkowski, de Sitter or anti de Sitter spacetime, 
and valid stable perturbation theory around these vacua. We secondly need to satisfy all actual 
experimental constraints as for ordinary GR plus we need correct IR cosmological behaviour 
without the need for dark energy nor a cosmological constant. We thirdly want our theory 
to have the least number of degrees of freedom possible and to be naturally connect to GR 
theory^. For example, Brans-Dicke theory [5] clearly passes the first and third tests whereas 
solar system constraints are rather restrictive [2]. 

1.2 Lovelock's theory 

In these notes we will restrict our attention to a metric modification of gravity that generalises 
GR in higher dimensions. Remaining tangential to GR (principles) we consider a theory 
£. = jC{M,g,V), whose field variable is a single symmetric metric tensor g endowed with 
a Levi-Civita connection V. We ask for a divergence free geometric operator on the right 
hand side of (2), since we know that matter obeys the conservation equation V^Ti^i, = 0. 
Furthermore, in order to bypass perturbative stability constraints for the graviton, we ask for 
second order field equations. These two properties are quite natural for our theory if we want 
to extend GR at the classical level but, we emphasize, not necessary at ultra-violet scales. 
Although higher derivatives generically introduce ghost degrees of freedom[7] ^ around the 
vacuum, [15] one may argue that these may dissappear having correctly summed the infinite 
number of higher order corrections. This is precicely the case in string theory which although 
is a ghost free theory of 2-dimensional surfaces embedded in 10 dimensions, at the effective 
action level, acquires (unphysical) ghost degrees of freedom because of the effective cut-off we 
impose. They arc in general cured by arranging for the appearence of the relevant Lovelock 
term [9] to the relevant order. 

In = 4 the only two derivative metric modification to Einstein's theory is the addition 
of a cosmological constant term! In other words, any higher order curvature invariant cither 
gives a pure divergence term, not contributing to the field equations, or, adds higher order 
derivatives to the field equations. In higher dimensions this no-go extension theorem to GR 

^Astronomical units are interesting since most tests of General relativity are at distance scales of the solar 
system. Hence extrapolation to lO'^^ scales bigger of such experiments can be sometimes unjustified or at least 
questionable. 

^The first and third requirements are not absolute but one needs to be aware at least when a theory does 
not validate one of these. 

''An exception to this rule are f{R) theories [8] since they involve only functionals of the Ricci scalar. 
These theories have been known since a long time to be conformally equivalent to scalar-tensor theories see 
for example [61] 
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is no longer true. It was the object of Lovelock's theorem [10], (see [11] for the D = 5 case) to 
prove back in the seventies that there exist theories containing, precise higher order curvature 
invariants, that actually modify Einstein's field equations (2) while satisfying V^T^i, = 0, 
in the face of modified Bianchi identities, and while keeping the order of the field equations 
down to second order in derivatives. The theory in question will be the subject of this 
brief study and gives in D = 4 precicely GR with a cosmological constant and in 5 and 6 
dimensions reduces to Einstein-Gauss-Bonnet theory (EGB). Lovelock theory in a nutshell is 
the generalisation of general relativity in higher dimensions while keeping the full generality 
of GR in D = 4. 

Following Lovelock's proof of the uniqueness theorem [10] (see also the neat derivation 
of [12] using differential forms) significant interest developped in these higher dimensional 
relativity theories in the eighties with motivations originating from string theory [13] and 
others originating from Kaluza-Klein cosmology [14]. Initial interest in string theory was 
trigerred by Zwiebach [15] who noted that second order corrections to the Einstein-Hilbert 
action, other than the Gauss-Bonnet invariant, introduced a graviton ghost when considering 
perturbations around flat spacetime. Effective action calculations of certain string theories [9] 
found that the leading (tree-level in gs) a' string tension corrections could give rise, modulo 
field redifinitions to this order, to the Gauss-Bonnet invariant. Several nice papers appeared 
uncovering properties and analysing exact solutions of EGB [16], [17], [18] while slightly later 
tackling full Lovelock theory exact solutions [19]. More recently there have been a few exact 
solutions discussed [20], [23] and some solution generating techniques developped (see for 
example [22] ) . Discussions on issues of energy, stability and the hamiltonian formalism have 
been carried out in [21], [26] [27]. 

Interest in Lovelock theory and in particular its 5 and 6 dimensional version, EGB theory, 
has attracted quite a lot of attention recently in the context of braneworlds (see Ruth Gre- 
gory's lecture notes [28]). Indeed from the braneworld point of view it would seem important 
to consider the general bulk theory rather than just GR in 5 or 6 dimensions and investigate 
if the 4 dimensional braneworld picture remained GR like. In a nutshell (we will uncover 
the details later on) the Gauss-Bonnet term in the bulk action is similar in nature as is the 
induced gravity term [30], [31] to be added to the brane action. Loosely speaking, it thus 
enhances GR type effects on the brane adding also quite naturally a UV modification to the 
usual one identified by the 5th dimension. Perturbation theory in the bulk is exactly the 
same around a maximally symmetric spacetime and the main difference are the boundary 
conditions on the brane which become mixed [29], similarily for those for induced gravity. 
In order to evaluate the correct boundary conditions, that give the braneworld gravitational 
spectrum, and hence determine the 4 dimensional gravity and stability of the setup, the im- 
portant difficulty one has to face, is finding the extension of the Israel junction conditions 
[32] in the context of EGB theory. In fact the junction conditions can be calculated directly, 
for each solution in question, by a careful calculation of the distributional tcrms^ [34], [35], 
[36] and [37] in the context of braneworld cosmology (see also [41]). The full covariant so- 
lution to the problem was first found by Davis, Gravannis and Willisson [38] where it was 
realised that careful variation of the bulk metric with respect to the boundary term to this 
theory, discovered by Myers back in the eighties [43], would give rise to the correct junction 
conditions. This is exactly similar to what happens when one considers careful variation of 

®As it was pointed out first in [33] distribution theory does not allow for ordinary multiplication and this 
can lead to erronious junction conditions. 
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the Gibbons-Hawking boundary term thus obtaining Israel's junction conditions. Using the 
junction conditions it was found that negative tension brancs induced tachyonnic instabihties 
to braneworlds [29] as well as important changes in the tensor perturbation amplitudes for 
braneworld inflation [39]. Braneworld cosmology was further studied with particular focus 
in inflation [40] and IR modifications [42], [27]. For codimension 2 braneworlds the relevant 
matching conditions were shown to give [44] precicely induced gravity terms on the brane 
plus extrinsic curvature corrections. However, up to now exact solutions or braneworld cos- 
mology have not been found. We will discuss briefly here how one can obtain the maximally 
symmetric braneworld solutions in the context of EGB [46]. The full matching conditions 
of Lovelock theory, irrespective of codimension were given in covariant formalism in [45] and 
recentl maximally symmetric braneworld examples to codimension 4 were found by Zegers 
[47] 

In this review we will therefore study the basic properties and important characteristics 
of Lovelock theory. In the next section we will begin by introducing the theory in differential 
form language as it is the most adequate way to recognise its nice features and why it has 
unique properties. After this geometric parenthesis, we will study important exact solutions 
of this theory concentrating on the case of static black holes and solitons. We will see that a 
generalised version of Birkhoff 's theorem holds as for OR, and we will then analyse the static 
black hole solutions, their thermodynamics and the solitonic solutions. After this we will 
discuss matching conditions of Lovelock theories and we will see that in this aspect Lovelock's 
theory is in principle a far richer extension to Einstein's theory in higher dimensions. Having 
done this we will discuss braneworld applications in codimension 1 and codimension 2. We 
will close off by looking at the Kaluza-Klein reduction of these theories and the type of 
scalar- vector-tensor theories they predict [48] . 

1.3 Basic definitions for Lovelock theory-differential form language 

Our aim in this section is to construct the higher order curvature densities which will be 
the building blocks of Lovelock theory and explain what makes them special. Indeed, in 
component language these will turn out to be precise but seemingly ad-hoc linear combinations 
or powers of the Riemann, Ricci tensor and Ricci scalar. We will thus use differential form 
language where we will see that they are indeed powers of the curvature 2 form with a precise 
and clear geometric interpretation (see also [69], [70]). 

Let {Ai,g,V) denote a D-dimensional spacetime manifold Ai endowed with a smooth ^ 
spacetime metric g. The connection V is taken to be a Levi-Civita connection. To every 
point P of spacetime M, we associate a local orthonormal basis^ of the tangent space TpM, 
{ca), with A = 1, D such that 

9{eA, en) = rjAB (4) 

Equivelantly the metric can be expressed as g = rjABG"^ ®6^= gabdx"'dx^ where the 1-forms 
Oa are precicely dual to the basis vectors e^, since ^"^(es) = 8^. The metric components in 
a coordinate frame dx"" (we use small case latin letters for a coordinate frame and upper-case 
latin letters for an orthonormal frame) are thus gab = O^O^rjAB where O^dx"" = 6^. The 
dual 1-forms 9^ form a natural basis of the vector space of 1-forms Q^^\TM). In turn the 

^By smooth wc mean hero metrics of at least regularity. This will be relaxed to piecewise when we 
look at braneworlds, allowing for distributional matter sources. 

''For a more complete account on the geometrical notions used here and precise examples see [77] 
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antisymmetric product of 1-forms 6 can be used in order to construct a basis of the higher 
order forms acting on TM. For any A;-form w in J^(fc)(TM), where {0 < k < D), can be 
written as, 

w = WA^-A,0^' A--- Ae^'^ (5) 
with WAi---Ak some smooth function. Following this simple recipe we can define a (D-k)-form, 

^M...A, = J^^eA,...A,A,^,...Aj^'^^ A---A9^- (6) 

where eyij Ao is totally antisymmetric in its D indices and ei2 - D = 1- This quantity is 

called the Hodge dual of the basis 0"^^ A • • • A 9'^'' of ^^'^^ (TM). It defines a dual basis of forms 
m ni^-''){TM). We can therefore write the Hodge dual of any k-form as, 

u = UA,...A,e^' A-'-AO^'- ^ *oj = u;^'-^''9X...A, (7) 

The wedge product of any form with its dual is a D-form, which is by construction proportional 
to the volume element of spacctime, which we note as 9*. Obviously for A; > -D all forms are 
identically zero. A useful identity is 

0^ /\ 0Ai...Ak = ^AA...Ak-i - ^Ak-i^Ai...Ak-2Ak + " " " + (-l)''~^'^Ai^A2...Afc • (8) 

Having constructed this tower of A;-forms, for < fc < we now need define two quantities: 
the connection 1-form and the curvature 2-form. The connection 1-form of M, lo^b which 
replaces the usual Christophel symbols in coordinate language, is defined by 

d9^ = -u^B A 9^ (9) 

since we have assumed a torsion-less connection. On the other hand, the spacetime curvature 
2-form is linked to the connection 1-form via the (second Cartan structure) equation, 

n^B = dw^B + u;^c A U^B ■ (10) 
The ambient curvature 2-form is related to the Riemann tensor components by 

n^B = \R^BCD9''A9'', (11) 

with respect to the spacetime Riemann tensor^. 

Langrangian densities appear under spacetime integrals; they are therefore by definition 
going to be Z?-forms. Hence in order to build such densities out of the curvature 2-form (11) 
and its powers, symbolically TZ'^, we need to construct i^-forms by correctly completing with 
the relevant Hodge duals (7). The higher the dimension of M the higher the possible powers 
of curvature we can consider. In differential form language it is straightforward to build this 
way the k-th Lovelock lagrangian density >C(jt) which is a D-form defined by, 

k 

jC-ik) = T^^'""' A • • • A 7^^'=^'= A 9Xs,...A,B, = A ^ ^MB,...A,B, (12) 

i=l 

*As a word of caution, the Riemann tensor components are given here with respect to the orthonormal 
basis and may differ, as is the case for stationary spacetimes for example, to coordinate basis components of 
the same tensor. 
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and we stress that k stands for the power of curvature. Clearly, £(o) tti^ volume element, 
giving rise to a cosmological constant whereas it is easy to check using (8) that 

£(,)=7^^lfilA0l^B^ = i^r, (13) 

is the Ricci scalar density and 

£(2) = n^'""' A 7^^2^2 A eXB,A,B, = {R^^'^'^Rabcd - ^r^^'Rab + R^) o* (m) 

is the Gauss-Bonnet density which we will denote by G. Note that for k > D/2, (12) vanishes 
so that if D = 4 say, then £(0), >C(i) and £(2) are the only terms present in the action (although 
>C(2); 'S'S wc will now see, turns out to be trivial). According to Lovelock's theorem [10], given a 
metric theory, C^^) the unique densities, made out of (A^, V,^) as defined in the beginning 
of this section, which allow for energy conservation and second order field equations. 



1.4 Lovelock densities and their geometric interpretation 

So what is special about these densities? The answer lies in differential geometry. Indeed, if 
D is even, for A; = D/2, the Hodge dual is trivial and the Lovelock density jO.d/2 reduces to, 

{D/2) 

= A T^^"^" ^AiBi...Ad/2Bo/2 ' (15) 
i=l 

This D-form can be recognised as the generalised Euler density for an even dimensional 
compact manifold Ai [78]. It is a geometric quantity whose integral over is a topological 
invariant (see for example [76]) 

^[^J= (4vr)^/W)! L'^^^/^^- ^''^ 

This relation has important consequences since it yields a relation between a geometric quan- 
tity involving curvature with the topology of the manifold M. This relation is familiar for 
the case of D = 2 where the A; = 1 density, ie the Ricci scalar or Gauss curvature, gives the 
usual Euler characteristic for 2 dimensional compact surfaces with no boundary. This famil- 
iar formula then relates surface geometry with a topological number defining topologically 
equivalent classes of surfaces, 

where h denotes the number of handles of M (see figure 1). The above relation gives in 
essence the Gauss-Bonnet theorem. This classification is also familiar from string theory 
where % is related to the string coupling gg = e^^ giving the string surface diagrams (rather 
than Feynmann diagrams as for point particles) for orientable strings. Variation of this 
quantity with respect to the local frame 6^ and use of the Bianchi identities on the curvature 
2-form, gives us that C.{pi2) is locally an exact form. This in turn tells us that variation of 
(15) gives no contribution to the field equations (supposing M. compact). In turn for D = 4, 
(l{2) stands for the Gauss-Bonnet density whose integral is now in turn the 4 dimensional 
Euler characteristic. 
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Figure 1: Any 2 dimensional compact surface can be continuisly deformed to one of its Euler 
classes parametrised by %. 



Therefore we see that the special feature of Lovelock theory is that its Langrangian densi- 
ties are dimensional continuations of the Chern-Euler forms^ which at lower even dimension 
are topological invariants. For completeness we note that for an even dimensional manifold 
with a boundary, the Chern-Gauss-Bonnet relation, is corrected by a Chern boundary form 
[78]. 

1 



X[M] 



(47r)^/2(D/2)! 



M 



<D/2) 



dM 



(D/2-1)) 



where the {D — l)-form $(d/2-i) reads [76, 78, 77], 



(D/2-1) 



D/2-1 
m=0 



D ■ {D -2)---2{m + l) 
l-3---{D-2m-l) 



Kl=Q 



2i-lM2i 



A 



N 



(18) 



(19) 



where K,^^ jv is a 1-form associated to the extrinsic curvature of the boundary dM. [69] , [45] . 
It turns out that the relevant boundary terms for Lovelock theory arc again the dimensional 
continuations of the Chern boundary forms which are multiples of the induced Riemann tensor 
on the brane and extrinsic curvature of the boundary. This is a rather useful statement when 
treating braneworlds since variation with respect to the local frame gives the codimension 1 
junction conditions [45] as we will see later on. The dimensionally continued boundary term 
associated to the Euler characteristic (13), /c = 1, is the Gibbons-Hawking boundary term [55] 
whereas the boundary term associated to the Gauss-Bonnet invariant is the Myers boundary 
term [43]. For details see [45]. 

The integral of the sum, for k < £>/2, of all Lovelock densities >C(fe) is the most general 
classical action for yielding up to second order field equations for the metric tensor. This 
is the Lovelock action. 



[(D-l)/2] 



where the brackets stand for the integer part. The first three terms of this sum 



Sd 



JM 



+ aiTZ^^ A e*AB + a2Tl^^ A 7^' 



^ ^ABCD H ) 



(20) 



(21) 



are respectively the cosmological constant, Einstein-Hilbert and Gauss-Bonnet terms, yielding 
the generalisation of the Einstein-Hilbert action in D = 5 and 6 dimensions. 



^By dimensional continuation we mean that we take k powers of the curvature 2-form and conveniently 
"multiply" by the relevant Hodge dual as in (12). 
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A variation of the action (20) including matter, with respect to the frame, gives the 
Lovelock equations 

[{D-l)/2] 

Yl c,k£ik)A = -2TABe*^ , (22) 

fe=0 

where is the k-th Lovelock (D-l)-form, 

k 

£{k)A = A ^""^""^ ^ 0AA,B,...A,B, (23) 
i=l 

and we have chosen the normalisation according to 

so that the equations of motion read, in component formalism. Gab + 
we get, 

where applying (8) and (11) we obtain in component language, 

Hab = ^{RefcdR^^"^ - ^RcdR"" + R') - 2RRab + ^RacR% + ^RcdR'ab " 2RaalxRb^^ 

the order 2 Lovelock tensor. In what follows we will explicitely set ao = — 2A, ai = 1 (expect 
in section where we will set a\ = Q and a2 = a for the Gauss-Bonnet coupling. 



• • • = Tab- For /c = 2 
(24) 



2 Exact solutions 
2.1 A staticity theorem 

We argued in the previous section that Lovelock theory is the natural generalisation of GR 
in higher dimensions. Which of the classical properties of GR remain true when we switch 
on the extra Lovelock densities. One main difference is that for /c > 2, ie once we allow for 
the Gauss-Bonnet term, (14), the Weyl tensor appears in the field equations. Therefore Ricci 
flat solutions are no longer solutions of Lovelock gravity if they are not conformally flat. This 
means in particular that construction methods quite common in higher dimensional GR are 
not going to give generically solutions for Lovelock theory. For example, it is not clear on 
how one can obtain a simple solution such as the black string [88], for Lovelock theory [22] 
and only a linear correction is known [95]. On the other hand we can question the status of 
classical GR theorems in Lovelock theory such as that of Birkhoff 's theorem. Is a version of 
this theorem still true in Lovelock theory? 
Consider the D-dimensional EGB action, 

R-2k + aG 

(25) 

where M is the fundamental mass scale of the D-dimcnsional theory, and G is the Gauss- 
Bonnet density (14). We note by 2A = —k'^{D — l){D — 2) the bulk bare negative cosmological 
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constant (2A = 0^(0 — 1){D — 2) is the positive cosmological constant), a is the Gauss-Bonnet 
couphng which has dimensions of length squared. We wih set a = {D—3){D — 4:)a to somewhat 
simphfy notation. This action is that of Lovelock theory in component language for D = 5 or 
6 dimensions. Let us for simplicity and without loss of generality [51] stick to D = 5 for the 
rest of this section in order to expose the staticity theorem. We will comment on the general 
result also involving an electric/magnetic field at the beginning of the next section. 

Suppose spacetime has constant 3-dimensional spatial curvature. This is the basic hy- 
pothesis we will make here. This is a slight generalisation from the common assumption of 
spherical symmetry in Birkhoff 's theorem. A general metric anzatz for the given symmetry 
is, ^ 

ds'^ = e^''^^''^B{t,z)-^{-dt^ + dz^) + B{t,z)T^ i -x'^^ 2 +X^dn]j_^ (26) 

where B{t, z) and u{t, z) are the unknown component fields of the metric and k = 0, ±1 is the 
normalised curvature of the 3-dimensional homogeneous and isotropic surfaces. We choose 

to use the conformal gauge in order to take full advantage of the 2-dimensional conformal 
transformations in the t — z plane. The field equations we are seeking to solve are found by 
varying the above action (25) with respect to the background metric or using (22) from the 
previous section and read 

£ab = Gab + ^gab-aHab = (27) 

where Hat is given by (25). 

It is rather useful to review and compare the equivalent system [71], [28] in Einstein gravity 
for q; = where, 

_ {D-S)A 

^ab ^j-^ 2^ y'^^' 



Pass to light-cone coordinates. 



t — z t + z , , 

u = ^, v = ^. (28) 



and take the combination R^t + R^z i 2Rtz = 0, one obtains the equations which read 

B,uu-2u^B,u = 0, (29) 

B^yy-2U^yB^y = 0. (30) 

Note then that these are ordinary differential equations with respect to u and v. They are 
directly integrable giving 

B = B{U + V) e^" = B'U'V' (31) 

where U = U{u) and V = V{v) are arbitrary functions of u and v, and a prime stands 

for the total derivative of the function with respect to its unique variable. We refer to (29) 
as the integrability conditions. Using a 2 dimensional conformal transformation which is a 
symmetry of (26), 

2 ' 2 

gives that the solution is locally static B = B{z) and the equivalent of a generalised Birkhofi"'s 
theorem is therefore true. Starting from a general time and space dependant metric, spacetime 
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has been shown to be locally static or equivalently that there exists a locally timclikc Killing 
vector field (here ^). By use of the remaining field equations we can then find the form of B, 
leading after coordinate transformation to a static black hole solution of horizon curvature k. 

Let us now take a ^ 0. In analogy to the previous case let us take the combination, 
Stt+£zz^^£tz = 0- On passing to light cone coordinates (28) we get after some manipulations 

(dB^/h^" + ISaKB^/^e^'' + 2aS „S ^) (fi - 2uuB,u) = 

(QB'^/^e^" + ISaKB^/^e^" + 2a5„B,„) (S,„„ - 2u^B^^) = (32) 

Note how the Gauss-Bonnet terms factorise nicely leaving the integrability equations (29) we 
had in the absence of a. 

The degenerate case where either B^u = or B^y = corresponds to flat solutions. For 
B^u / and B^^ / the situation is clear: either we have static solutions and the staticity 
theorem holds as in the case above or we will have 

^2. ^ MBl - Bl) 

952/3(52/3 + 4aK) ^ ' 

Let us briefly examine the latter case, that we will call Class I solution [17], [37]. The 
two remaining field equations E-^-^ = and £tt — £zz = are solvable iff we have the simple 
algebraic relation, 

Aak^ = 1 (34) 

This is quite remarkable: if the coupling constants obey (34) then the B field is an arbi- 
trary function of space and time; in other words the field equations do not determine the 
metric functions. Therefore strictly speaking Birkhoff 's theorem does not hold for non zero 
cosmological constant. Setting B{t,z) = R^{t,z) the class I metric reads. 



ds' = '\ ^\ -dt' + dz') + R'[ rr^^ + x'mi ) (35) 

2a 



This solution has generically a curvature singularity for = ±-R,t. The Class I static 
solutions are given by, 



.2 _ A{Rr ..2 , dR' , „2 f dx' 



ds^ = -^-^dt^ + — ^ + R^ + x'dnjj (36) 

'^+2^^ '^'^2^ V-^'^A / 

with A = A(R) now an arbitrary function of R. 

In order to obtain t and z dependent solutions it suffices to take the functional R to be 
a non-harmonic function. Take for instance R = exp(f(t) + g{z)), with / and g arbitrary 
functions of a timelike and spacelike coordinate respectively. Let us also assume At = for 
simplicity, the Class I metric in proper time reads. 



ds' 



l + 2a/2 Vl-^X^ A jyy 



^°Note however that for a non-zero charge Q and spherical symmetry (ft = 1) Birkhoff's theorem is always 
true as was first shown by Wiltshire [18] 
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Note here again that / is an arbitrary function of time and g an arbitrary function of space. 
The fine tuning relation between a and k actually corresponds to a case of enhanced symmetry 
often refered to as Chern-Simons gravity (for odd dimensional spacetimes [58]). 

On the other hand if (34) does not hold then Birkhoff's theorem remains true in the 
presence of the Gauss-Bonnet terms i.e. the general solution assuming the presence of a 
cosmological constant in the hulk and 3 dimensional constant curvature surfaces is static if 
and only if (34) is not satisfied. In this case the remaining two equations give the same 
ordinary differential equation for B{U + V) which after one integration reads, 

(B' \ ^ 

+ =%ix (38) 

where ^ is an arbitrary integration constant. Then by making B the spatial coordinate and 
setting r we get the solution discovered and discussed in detail by Boulware-Deser 

[16] {k = 1) and Cai [52] (k = 0, -1)," 

ds- = -V{r)dt- + ^ + + X'd^]) (39) 

where V{r) = k + ^[Irt ^1 — 4aP -|- 4^], and ji is an integration parameter related to the 
gravitational mass. The maximally symmetric solutions are obtained by setting /x = 0. We 
will analyse in detail these solutions in the following section. 

To close off notice how (34) is a particular 'end' point for (39) since the maximally sym- 
metric solution is defined only for 1 > 4q;A;^ (for a < there is no such restriction). We 
can deduce in all generality that for 1 > Aak'^ there is a unique static solution (39). When 
(34) is satisfied and fi = 0, the two branches coincide (we have an infinity of solutions) and 
V = K i: \J nja + ^, is then a particular Class I solution (34) very similar to the BTZ 
3-dimensional black holes. For 1 < 4aA;^ no solutions exist. We will come back to the 6 
dimensional version of these black holes in a moment. Before doing so let's see briefly how 
they generalise for the full Lovelock theory. 



2.2 Lovelock black holes 

The staticity theorem we evoked in the previous subsection is generalised without major 
difficulty for the general Lovelock theory in arbitrary dimension and in the presence of an 
Abelian gauge field [51]. That is for the theory involving a Lovelock action (20) with a 
Maxwell gauge field, 




we obtain from the field equations that apart from the pathological cases of Class 1 there is 
a unique solution given by [19] 

ds- = -Vir)dt^ + ^^+ r\^, + x'dnU) (41) 
^'^We have kept the same label as in (26) for the rescaled time coordinate . 
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where the electric field strength is F = ^^^2(0-2) dt A dr. The metric potential reads V = 
K — r'^f{r) where / is a solution of the fc-th order algebraic equation, 

P(/) = E"'/' = ^-;:& (42) 

1=0 

where k = [^y^] is the order of the Lovelock theory. It is clear that the higher the dimension 
of spacetime, the more the terms in the Lovelock action, and hence the higher the order of 
the equation (42). The k possible roots of the polynomial P{f ) for k = 1 actually give us 
the maximally symmetric vacua of the theory [19] . Note that even in the absence of a bare 
cosmological constant these vacua can be flat or of positive or negative curvature and their 
magnitude depends on the normalised Lovelock coupling constants, 

"0 1 . _ 

a^{D-l){D-2y 
21 

— WiD-n),ioT 1>1 (43) 

Positive roots correspond to dc Sitter vacua whereas negative roots to anti de Sitter vacua. It 
is interesting to note that in the presence of a bare cosmological constant ckq 7^ maximally 
symmetric vacua may not exist at all. For zero bare cosmological constant however the flat 
vacuum is always solution. Solutions of (42) have complex horizon structures and have been 
analysed by Myers and Simon [19] for Q = 0. 



ao = 
ai = 



2.3 Einstein-Gauss-Bonnet black holes 

For simplicity, let us now truncate Lovelock theory at = 2 ie neglect higher order terms 
other than the Gauss-Bonnet invariant. We will follow for most of this section the analysis 
of [19]. We therefore consider the action (40) for up to A; = 2, 

Si^ = 5(2) -\l d^xV^Fa,F^' (44) 

Take a {D — 2) dimensional space of maximal symmetry and therefore of constant curvature 
parametrised by k = 0, — 1, 1, 

= h^dx^dx'' = + x'dnl^_^) (45) 

These are the maximal symmetry spaces we considered in the previous section and express 
the constant curvature geometry, normalised to k, of the horizon surface. If on the other 
hand we perform a careful Wick rotation to h'^ we can construct h^, which is of Lorentzian 
signature, and gives spacetimes which are {D — 2)-dimensional sections of Minkowski, adS and 
dS. Therefore, the staticity theorem of the previous section tells us that any D dimensional 
spacetime metric admitting D — 2 dimensional maximal sub-spaces (or sub-spacctimes h^) 
which is solution of the field equations emanating from (25) is locally isometric to, 

ds^ = -V{r)dt^ + — - + r'^hl^dx^dx" (46) 
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with electric field strength 



47rr2(^-2) 

or, modulo a double Wick rotation, to 



«2 

^ dtAdr (47) 



ds^ = y(r)d02 + + r^hj^^dx^'dx" (48) 
with magnetic field strength 

^ = i;^;:fc5)^^^^- (49) 

In the latter case the theorem gives a local axial Killing vector do and concerns locally axially 
symmetric solutions. This case has been explicitely treated for the case of general relativity 
= 1 in [59]. Wc will come back to this simple yet powerful result to give the maximally sym- 
metric cosmic string metrics in D = 4 and maximally symmetric braneworlds of codimension 
2 in Z) = 6. 

Since we want to search for black hole criteria wc concentrate on the static case for the 
rest of this subsection. The potential reads (as can be easily verified from (42)), 



V{r) = f 

with integration constants, 



2a 



l + 6yi + 4a(a2 + -A^--Ji^) 



(50) 



^2_ 167rGM 

^ 2TriD-2){D-3)' ^~{D-2)^^ ^ > 

where q is the charge, M is the AD or ADM mass of the solution and is the horizon area. 

First notable fact is the the ambiguity of the vacuum which is parametrised by e = ±1 and 
gives rise to 2 branches of solutions. Indeed setting = = 0, k = 1 gives us the possible 
vacua of the theory (25). In fact wc note that if wc set the bare cosmological constant to 
be zero a = we note that we do not only obtain the fiat vacuum. For e = 1, we actually 
asymptote anti de Sitter space for a > and dS space for a < 0. We will refer to this vacuum, 
as the Gauss-Bonnet branch. Indeed the effective cosmological constant is Ag// = 
and hence a plays this role without a bare cosmological constant term in the action. Unlike 
what was first argued in [16] this branch is not, at least, classically unstable [49] in the sense 
that one needs to add positive energy to the system for it to roll off to a positive mass black 
hole solution. We will here caution the reader that this branch is still dangerous for stability 
and a further careful analysis is still needed to answer the question of the physical relavance of 
this branch [50]. It is an intriguing fact however, that one can have an effective cosmological 
constant from higher order curvature terms as used in [60]. Any solution of this branch will 
not have an Einstein theory limit, although there maybe a relevant de Sitter or anti de Sitter 
Einstein type solution mimicking (46) as wc will see in a moment. The e = — 1 branch gives 
the usual Minkowski vacuum in the absence of a bare cosmological constant A and for large 
r the solutions resemble the asymptotically Einstein black holes with the relevant mass and 
charge parameters. This branch therefore is perturbatively connected to Einstein theory and 
we can consider the action (44) as an effective action including a higher order correction as 
in actions for closed strings [9]. For the Gauss-Bonnet branch e = 1 notice that this is no 
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longer the case and small a yields an enormous effective cosmological constant. Hence it is 
in this case an erronius statement to consider (44) as an effective action and one needs exact 
solutions rather than perturbative ones. This fact makes this branch either totally irrelevant 
or far more interesting than the Einstein branch since this is where we can expect novel effects. 
Given therefore the fact that the e = — 1 case closely follows known solutions we will check 
out mostly the Gauss-Bonnet branch for novel effects. 

The "Chern-Simons" limit is obtained for 4aA;^ = —Aao? = 1 (for a positive and negative 
respectively). The two branches then become one and at this limit the combined Gauss- 
Bonnet and cosmological constant are of the same order as the Einstein term on the right 
hand side (whose coupling we have normalised to 1). For zero charge the potential simplifies 
to the simple function 

2 

V TT) 

V(r) = i^+t: 7— (52) 

where m = —^J^ is the mass integration constant. For D = h the potential is quite similar 
to the BTZ black holes of 3 dimensions since the mass plays the same role as curvature of the 
horizon k. Furthermore note that we have a weaker gravitational force than in the Einstein 
case. Indeed set D = Q and a planar black hole for a > exists with > (2ma)^/^ whereas 
the Einstein horizon for the same mass m is its square root. Hence for the same mass the 
Chern-Simons black hole is of squared horizon radius compared to the Einstein one. Seemingly 
the closer we approach the Chern-Simons branch the milder is the curvature singularity in 
the bulk spacctimc an interesting fact to keep in mind. 

So much for the branches of solutions. Let us now check out the curvature singularities 
and horizons. Generically, there are two possible singularities in the curvature tensor, the 
usual r = 0, but also a branch singularity at the maximal possible zero of the square root, 
say r = ri . We have that r = ri is solution of 

^2(D-2)(i ^ 4^^2) ^ 4c,(^ri?-3 -Q^)=o (53) 

and whenever n > this is the singular end of spacetime (46). We have a black hole solution 
if and only if there exists r spacelike with r = such that V{rh) = and Vh > rmim where 
"""min = max{0,ri). Indeed the usual Kruskal extension 

dv 

dv± = dt± — - (54) 
V{r) 

gives that (w+,r) and (f-,r) constitute a regular chart across the future and past horizons 
of (46) as in GR. It is straightforward to show the following criterion : r = r/j is an event 
horizon for r > spacelike iff, 

• fh ^ ffnin 

• €{2aK + rl)<0 

• Vh is root^^ oi pa{x) = (— a^a;^(^~^) + kx'^^^~^^ + an^x'^^^'^^ — iJ,x^~^ + Q'^)sign{a) 

Notice therefore that for e = 1 we have r| < —2aK ie the event horizon is bounded from 
above. This immediatelly yields that for «; = there are no black hole solutions in this 



^For r spacelike we need p{x) > away from the horizon for a > and the contrary for a < 
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branch! Secondly, note that Pa=o is the usual polynomial for A-RN black hole of Einstein 
theory in D dimensions. In particular, a couples only to the horizon curvature. Hence for 
e = — 1 and «; = (the Einstein branch) the horizon positions are the same as for the GR 
planar black holes. 

Let us now focus on some particular solutions for the Gauss-Bonnet branch (again we 
advise the interested reader to see the nice analysis of [19]). Thus take e = 1, a = 0, k = 1. 
We have < —2a, hence take |a| = —a > 0. We also take < in order for (53) to be 
strictly positive and to have a correct definition of mass. Indeed the solution asymptotes 
a de Sitter Schwarzschild solution with positive AD mass for negative /x. Now given that 

Pa{x) = —x^{x^ — x\a\ — fi) we find Xmin = have two event horizons as long 

as < fi < 0. The bigger the mass of the black hole the further we are allowed to 

stretch the a parameter of the action (see figure). The solution has the same structure as 
Schwarzschild de Sitter with an event horizon and a de Sitter horizon. If we switch on we 
can get a 3 horizon structure similar to the de-Sitter RN black hole. When the two horizons 
come together at Xmin we are in the extremal case of the Nariai solution. 

Consider now a hyperbolic horizon, in other words n = —1. Take Q = a = and e = 1. 
We have = — 4a/x and > —2a. Therefore for a > and /x > we have trivially the 
first two conditions verified. Furthermore, Vh is root of —x^ + ax — fi which is exactly the 
same polynomial as above for a — > —a and ^ -ji. Therefore in this case we will have a 
hyperbolic black hole with 2 event horizons. 



2.4 Thermodynamics and a geometrical explanation of the horizon area 
formula 

Our aim here is to evaluate the entropy of the black hole solutions descussed in the previous 
section. We therefore start by evaluating the mass and temperature. The mass of the black 
hole can be easily expressed in terms of the horizon radius, by. 

In order to calculate the tempcratTirc of the black hole we follow the standard prescription. In 
summary: we start by Wick rotating the time direction t ^ iO. The resulting curved manifold 
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of Riemannian signature has an axial Killing vector de at its origin situated at r = r/i as we 
will see in a moment. We then impose periodicity, say IT, in order for the angular coordinate 
to be well defined. As a result the Euclidean quantum field propagator, with the imposed 
periodic boundary conditions, describes a canonical ensemble of states in thermal equilibrium 
at a heat bath of temperature T = [62]. 
Indeed consider t ^ iO oi (46) . We have 

ds' = V{r)de' + y-r + r^hl^dxi'dx" (56) 



with (magnetic) field strength 



l + 6yi + 4a(a2 + -ii^ + -^: 



(57) 



P = ^;^^)d(^^dr (58) 

and the metric is then of Euclidean signature for r > Vh and p = iq (P = iQ) is the magnetic 
charge. Let us put the charge equal to zero for simplicity. The potential V{r) admits at least 
one and at most 2 zeros for < r_ < r_|- where r = r_|_ can be taken to be infinity for the 
single horizon case. Taking x'^ = constant and expanding around the zeros of the potential 
V we get, 

ds'^ ~ (iK±) P±dO^ + dpi (59) 

with radial isotropic (or cylindrical) coordinate p± = ^ ^('^^±) "_ Hence the periodicity condi- 
tion reads, 

47r Air 

e^e + — ^9 + — (60) 

Generically we will have conical singularitites. In fact the only "finite"' regular solution 
exists when, V^'_ = V^^ ie for a T = temperature black hole (Nariai limit) and it corre- 
sponds to a regular gravitational instanton. In fact, one has to consider a simple coordinate 
transformation to show that space is of non zero volume and corresponds to x 5^ check 
out [56], [46]. When we add charge it is a question of algebra to show that we can obtain 
warmer, T ^ instantons^^. Gravitational instantons describe the extremal spacelike path 
between two non-causaly connected regions of gravitational solutions [82], in other words 
their quantum tunneling. In the case here the instanton describes decay of the vacuum into 
a pair of accelerated black holes [46] as in usual GR (for instanton solutions using boundaries 
see [57]). The decay rate of course would have to be calculated here from scratch since it 
corresponds to a novel theory [46]. If it is enhanced it corresponds to a quantum instability 
of the vacuum. When r+ = oo space rounds up smoothly and regularity is only imposed at 

After this brief parenthesis let us now calculate the temperature of the black holes (46). 
The period is 11 = 4r^. It is easy to show the formula, 

V'y._|_ 



v =1- 



r^"^(2aK + 4) 



(61) 



^^I thank Renaud Parentani for pointing this out 
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The heat bath described by the thermal propagator is of temperature T = 11 ^. Therefore 
the temperature is found to be (using (55)), 

^ ^ {D - l)fcVg-^ + {D- 3)Kri-' + {D- 5)aK'ri-' 

Notice that the Gauss-Bonnet term a couples only to the horizon curvature and the tem- 
perature is the same for e = ±1. Also note that a Lovelock planar black hole will have 
same temperature as an Einstein one [52]. For the entropy we use here the standard recipe, 
dM = TdS (see also [26] lor a direct calculation yielding the same result) . Then in turn. 



J Jr ■ <)fh 

* 777, 7, 77, ' * 



- if:!^^ 2a^<£-2) 

" 4G {D-A)rl ' ^^^> 

In order to evade erronius conclusions it is important to note that rmin = max{ri,Q) [81], 
since in Lovelock black holes we can hit a singularity before reaching r = 0. Furthermore we 
see that the entropy is not equal to the area of the black hole horizon; we pick up a correction 
from the induced curvature of the horizon surface. One can understand this fact geometrically 
using the general formalism of Iyer and Wald, [80]. The entropy calculation in this case [81] 
gives the result, 

1 I — 

S = — dx^-^ V V(l + Sai?^"*^) (64) 



where the tilded terms correspond to geometrical quantities of the horizon surface, (r, 9) =constant. 
The result is exactly the same as the effective action obtained for a codimension 2 matching 
condition (85) as we will see in the forthcoming section^^.. This is not too surprising with 
hindsight: In both cases we impose similar regularity conditions which have to do with the 
temperature and hence the periodicity of the manifold. In one case conical singularities are 
accounted for by the presence of branes of given tension whereas for black holes this preci- 
cely gives the temperature of the horizon. The leading Lovelock correction appearing as an 
induced curvature term is therefore nothing but the extended Euler density of the horizon 
surface. On the other hand the Einstein term yields the tension of the horizon. Therefore 
the failure to obey in this case a horizon-area formula is most natural and a simple geometric 
consequence. The entropy of the Lovelock black holes follows the matching condition formula. 
For example if we were to take a D = 8 dimensional Lovelock black hole we would get in 
addition the 6 dimensional Gauss-Bonnet correction of the horizon surface in the entropy 
formula namely, 

' ^(1 + 2a^^"'^ + aaG*"*^) (65) 




We will come back to these formulae in the next section. 



2.5 The Lovelock solitons and the rrictximally symmetric cosmic strings 

So much for the Euclidean version of (46) and the thermodynamics (for details check [81]). 
We now Wick rotate the D — 2 horizon sections h^^ to the Lorentzian maximally symmetric 

^^Here the extrinsic curvature quantities appearing in (85) cancel because the horizon surface is maximally 
symmetric 
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sections hj^^ to construct soliton solutions (48) (the same obviously holds for the Lovelock 
case (41)). Since the soliton is of axial symmetry wc have r_ < r < r+ where r± are the 
possible zero's of V. One can always locally go to the cylindrical anzatz (59) with p± the local 
cylindrical coordinate. The soliton is of infinite proper distance p when r+ ^ oo. A very nice 
and simple interpretation of these solutions and a simple application of the statiticity theorem 
resides in D = 4 general relativity. In this case the Lorentzian 2-dimensional sections in (48) 
have precicely the geometry of a maximally symmetric cosmic string, of internal constant 
curvature given by k, situated at r = r± once we allow for some deficit angle in 9. These are 
the only solutions describing de-Sitter, fiat or anti de-Sitter cosmic strings in a cosmological 
constant and (possibly) charged background. We will not analyse these here but as an example 
it is easy to note that in the presence of a bare negative cosmological constant the straight 
cosmic string bends the ambient spacetime (this was first noted by Linet in [64]). This is in 
complete contrast to the fiat string which is embedded in a locally flat spacetime. This is 
because the resulting geometry is of non-trivial Weyl tensor, 

precicely denoted here by the presence of a non-trivial p parameter. In (66) (t, z) are the 
string coordinates, the angular deficit is given by ^ = 27r(l — 3) and the string is of linear 
tension T = ^^q^ ■ In other words a straight cosmic string mathematically at least corresponds 
to a double Wick rotated black hole solution. Another way to understand the above result is 
that in pure adS spacetime we cannot slice the geometry in a cylindrically symmetric anzatz. 
It is a simple exercise to classify the solutions, the case of de-Sitter also presenting some 
unique characteristics due to the compactencss of its spatial sections (cosmic strings will 
appear in pairs in this case). Furthermore, for k = 0, D arbitrary and negative cosmological 
constant we get the Lovelock ads soliton analysed in [79] for the case of Einstein gravity. In a 
forthcoming section we will see that in D = 6 dismensions the metrics (48) will describe (with 
the help of suitable matching conditions) the generic bulk geometry of max;imally symmetric 
4 dimensional braneworlds [46]. 



3 Matching conditions for distributional sources 

Most recent applications of Lovelock theory concern braneworld physics. In the braneworld 
picture our 4 dimensional universe is part of a higher dimensional manifold and whereas 
matter is strictly confined on the braneworld, gravity propagates, according to the equivelance 
principle, in all dimensions. We therefore want to describe, the motion or the evolution, of a 
self-gravitating submanifold which for simplicity we take to be infinitesimaly thin. In other 
words, we suppose that matter is confined on the braneworld via a Dirac distribution of 
dimension equal to the codimension N = D — 4 of the braneworld. Given that Lovelock 
theory is the general metric theory with second order field equations we can except some 
junction conditons quite similar to those of Einstein thcory^^. It turns out that in Lovelock 
theory we can even do better, at least in mathematical terms, and define matching conditions 
for higher codimension than that possible for Einstein theory [63]. 

^® Fourth-order field equations would mean that distributional constraints would be imposed on the continuity 
of certain third order directional derivatives which would mean a discontinuous limit to pure Einstein theory. 
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Consider a p-brane S embedded mD=p+l + N dimensions and suppose that E carries 
some localised energy-momentum tensor 

Tab=^ \ (67) 

The A^-dimcnsional Dirac distribution Sy. on S signifies that the brane is of zero thickness. 
The question we address is: What are the equations of motion for this self-gravitating p- 
brane sourced by the distributional energy momentum tensor (67)? For codimension N = 1, 
the answer is given by the well-known Israel junction conditions [32] where, if the induced 
metric is continuous, a discontinuity in the first derivatives of the metric accounts for the 
Dirac charge in (67). Israel's junction conditions describe adequately a wide variety of GR 
problems and only when we study junction problems of lesser symmetry do we run into 
shortcomings because we can no longer fullfil the continuity condition (for example we cannot 
match Kerr to flat spacetime) . If the codimension is strictly greater than 2 then there are no 
distributional matching conditions in Einstein gravity. Thus, a finite thickness braneworld 
is needed in order to obtain non-trivial self-gravitating equations of motion for S and the 
resulting equations of motion will generically depend on the regularisation scheme. This is 
hardly surprising. We know for example, that far away from a gravitational source such as the 
sun, and by virtue of Birkhoff 's theorem, we can approximate conveniently its gravitational 
field using a 3 dimensional Dirac distribution. What we mean by far away is precicely the 
Schwarzschild radius of the source in question which is about 3km and which is far smaller 
than the actual size of the sun. In Einstein theory even when the codimension is equal to 
2 (for D = A say), one only knows the self-gravitating field of a straight cosmic string, ie 
one induced by a pure tension matter tensor, which gives an overall conical deficit angle [63]. 
Hence even for codimension 2 distributional matching conditions only give a brute picture of 
the gravitational field and break down when the braneworld or defect is of lesser symmetry. 
The essential point for the discussion is the codimension of spacetime A'^, ie the number of 
spacelike or timelike, vectors defined normal to the brane (for lightlike junctions see [65]). If 
the codimension is 1 as it is for the usual junction conditions the p-brane is a hyper surface 
splitting the bulk in two and there is no geometry (other than that of real line) in the normal 
directions. Once > 1 we have non-trivial geometry in the normal sections. 
Lovelock's equations for the distributional energy momentum tensor (67) are 

[(D-l)/2] 

J2 ^k£(k)^. = -2Si,,e*''Sj: , (68) 

k=0 

Our task here involves finding the distributional part of the LHS to be matched to the induced 
energy-momentum tensor S^i,. We will refer to such geometric distributional terms as those 
carrying Dirac charge. It turns out that not all Lovelock terms can carry a Dirac charge, we 
already know this from Einstein theory which is Lovelock theory in D = 4. Indeed we will 
find the simple inequality, 

N/2<k< , (69) 

which selects in particular the Lovelock bulk terms >C(jt) that can carry a Dirac charge. Indeed 
for Z) = 4 we have N/2 < k < 1 which tells precicely that an Einstein theory, k = 1 can carry 



^® Greek letters run through brane coordinates while capital Latin from the begining of the alphabet through 
bulk coordinates 
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up to codimension 2 Dirac charge. If D > 4 and we allow for the higher order Lovelock terms 
we see that we can go to higher codimension. 

Without attempting to give a full proof of this result, given in [45], we can however 
understand it geometrically in the following way: seperate the geometry in normal and parallel 
sections as we commonly do in the Gauss-Codazzi formalism. Let e^^ be the p + 1 unit vectors 
that are everywhere tangent to the brane, n/ be the A'' unit vectors that are everywhere 
normal to E with the label M denoting the radial normal vector. Similarily, split locally the 
bulk 9^ into tangent 1-forms 9^^ and normal 1-forms 9^. One can then deduce the first and 
second fundamental forms of the brane, respectively; 

h = Vf^uO^^e" (70) 
Kif^u = g{'Ve^nj,e^) (71) 

describing the induced metric and extrinsic curvature of S. The Gauss-Codazzi equations 
can also be written in form formalism which is extremenly useful in this context. Indeed the 
parallel projection along the brane coordinates gives the Gauss equation. 



7^;/^ ^ ^0' = - IC'j/j A /Cj/, , (72) 



where 

^//m = ^^^'^ ' 

is the extrinsic curvature 1-form and Jl'f, is the induced curvature 2-form of the brane, 

//I' _ ... 
associated with the induced metric This geometrical identity relates the background 

curvature with respect to the induced and extrinsic curvature of S. 

We look for geometric terms in the //z^-components on the LHS of (68) that can carry a 

Dirac charge. The situation is inherently different for odd and even codimension. Indeed for 

hypersurfaccs (codimension 1), the distributional charge in Einstein theory is provided by a 

discontinuity of the extrinsic curvature of the metric Kj^^^. Thus we relate local geometry of 

the brane to matter. In even codimension defects, like cosmic strings (codimension 2), the 

normal sections to the brane have some non-trivial topology that can give rise to distributional 

terms. For the case of the straight cosmic string we have, 

ds"^ = -dt^ + dz^ + dp^ + L^{p)d9^ + dp^ (73) 

and Einstein field equations are solved for L{p) = p everywhere but at p = 0. We therefore 
set 

L'{p) = 1, p^Oand L'{p) = /?, p = 
The distributional part of (68) is then integrated over the normal section and reads 

/ / dist{L")rdr d9 = -8TrG4 —6^^^rdrd9 (74) 
Jo Jo Js 27r 

Given that dist(L") = [L']^^^) = we get the well-known result 27r(l-/3) = -SttGaT 

relating here topology (rather than geometry) to matter. It is important to note that the 
cuvature of the normal section is precicely given by L". 

In each case we define locally a normal section, T,±, at each point of the brane, using the 
congruence of the normal vectors with specific regularity properties. We then integrate the 
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juz^-equations of motion over an arbitrary such section Sj^ • Since the Dirac distribution is by 
definition independent of regular isation ie to the microscopic features of S, the only terms 
that can contribute are those independent of geometry deformations over T,± as we take the 
limit of zero thickness. This boils down for codimension 2 to working out the distributional 
part as in (74). Mathematically, the only terms in (68) having this property are proven to be 
locally exact forms on T,±, basically made out of powers of extrinsic curvature or curvature 
of the normal section as for the example of the cosmic string above^^. 

What multiplies these locally exact forms are now equated to the energy- momentum tensor 
of the brane and depend on the parallel geometry (since the normal geometry is integrated 
out) using essentially Gauss's equation (72). Since higher order Lovelock terms, k > 1, 
contain higher powers of the curvature tensor, we can expect the appearence of induced 
curvature and extrinsic curvature terms from (72). Indeed the relevant terms are sums in 
< j < min{[N/2], k-N+ [N/2]) of 

/k-N+[N/2]-j 

-(k.>= A ^rlM A ^?/V|a^;1,,, (75) 




which are powers of the projected Riemann curvature two form TZ// and the extrinsic curva- 
tures fC//. 

For the case of codimension 1 the junction conditions read, 

[(D-l)/2] 



- E 



Sakkl 

k=0 



a" 



2S/(P)C// . (76) 



and the Dirac charge will be provided by the jump of c^jy^.)^- Here 1 < A; < [{D — l)/2] 
and the higher the dimension the more the Lovelock charges that contribute. For example in 
p = 3 we get in turn for the k = 1 Einstein and k = 2 Gauss-Bonnet terms, 

= '^ic'//N^o;^/ = -2{K;-s%K)ef 

= -4{3j;:-J6;-2P\^K'n0f, (77) 



where we have set 



PtivXp = R/xuXp + '^Ru[\gp\ii - '^Rp{\gp\v + Rgti[\gp]u 

and we have dropped the label N for the extrinsic curvature components. Replacing into 
the matching conditions (76) and taking the corresponding left and right limits one gets the 
well-known junction conditions for Einstein [32] and Einstein-Gauss-Bonnet gravity [38]. The 
equations of motion can also be derived from the boundary action, after variation with respect 



^^In codimension 2 the above statement boils down to the famous Gauss-Bonnet theorem. In higher codi- 
mension the relevant forms are precicely the Chern-Simons forms, [78] 
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to the induced frame 9'^ on the left and right side of the brane. In differential form language 
the action is obtained trivially from (77) by removing the free index, 

5e = 2ai ^ /CJ/^ A e*/ + 4a2 ^ /CJ/^, A (nf - ^iCj^^ A /C)jv) A 9lf, (78) 

and agrees with the Gibbons-Hawking [55] and Myers [43] boundary terms. It is obvious now 

that for Z? = 7 a 5-branc will have an extra term contributing to the junction conditions which 
will be of 5th order in the extrinsic curvatures etc. This can be understood intuitively since 
higher order Lovelock densities involve higher powers of curvature. For the case of higher odd 
codimension see [66]. 

Let us now concentrate on even codimension N = 2n. We will assume that the parallel 
sections are regular. Then when integrated out over T,±, the normal sections quite naturally 
yield the charge or topological defect of (x — P)Area{S2n-i)- This is similar to the defect 
angle for the infinitesimal cosmic string in 4 dimensions (73) when we set x = 1; where x is 
the Euler characteristic of the normal section^^. In higher codimension one has to allow for 
non-trivial Euler character in T,± in order to obtain removable type singularities [47]. Due to 
their topological character, we call these topological matching conditions [66]. The matching 
conditions read. 



{n,k)n 

k=0 



ix - /3)Area(52„-i) ^ ^ - (P) = -2S;{P)9t// , (79) 



where we have set = 2^" ^fc!(n — l)\ak/{k — n)\ and the smooth parallel forms a^^^ are 




(80) 



Finally, k = k — n, which ranges between < k < [|] , will turn out to be the induced 
Lovelock rank of (80). The parallel forms (80) will dictate the dynamics of the brane. Their 
expressions involve powers of the projected Riemann tensor TZ// on the brane and even powers 
of the radial extrinsic curvature /Cjv of S. To see this, consider for each k the first term in 
the sum, j = 0, which reads, 

^// _ I A -T'^i'^i 1 A a*// 



(fe,0)// \ / \ // I Mli^l-Aj^^fc 
\l=l J 

Note the similarity to the bulk Lovelock densities, (23). Clearly, using (72), we see the 
appearance of induced Lovelock densities involving fi// accompanied by even powers of the 
extrinsic curvature K.N. The former describe the induced quantites of the brane and the latter 
how S is embedded in the bulk. In particular, for = we will have a pure tension term, for 
= 1 an Einstein term and for k = 2 a Gauss-Bonnet term with extrinsic curvature terms. 
Note that the highest rank on the brane originates from the highest rank Lovelock term in 



*I thank Robin Zegers for pointing this out 
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the bulk. This agrees with the fact that in Einstein gravity, k = 1, there are no matching 
conditions beyond n > 2. 

Let us compare our result with the work of [44]. Thus, let us consider the case of a 3-brane 
embedded in 6 dimensional spacetime, for which there are only two terms in (80), 

<o)M = <i)M = ^7/^<V (82) 



Thus, using (8), (11) and (72), we obtain from (79) 



27r (1 - P) {-ai V + 4a2G(::^'^^ - 4a2 W^,} = 



Sfiu , (83) 



where [44] 

K = KnK^^n - KInK^n - [Kl - k;^K^) . (84) 

Note that the bulk Einstein term, A; = 1, only allows for an effective cosmological constant on 
the brane, while the Gauss-Bonnet term, = 2, induces the Einstein tensor for the brane's 
equation of motion. This equation is similar to the one found in [44], the difference being 
that here, the extrinsic geometry is supposed perfectly regular. Indeed mathematically there 
is no reason to suppose that the extrinsic curvature has a jump, the topological defect car- 
rying the necessary charge in complete analogy to the case of the cosmic string. Aditionally, 
these matching conditions provide the maximal regularity for the bulk metric. Furthermore, 
since the induced and extrinsic geometry are smooth, the equations of motion can be seen to 
actually originate from a simple action taken over S as long as we suppose that (3 = constant. 
In other words, the degrees of freedom associated to the normal section are completely inte- 
grated out giving an exact action for the brane's motion. In differential form language, it is 
straightforward to read off the langrangian density in question. Literally, taking out the free 
index for the charge in (81) and using Gauss equation (72), we get 

^(p=3,n=l) ^ 2nil-P) J^(^di9*// +d2inj A 9*4))+ J^Crnatter 

= 27r(l -p)J^V^ (ai + a2(i?'"'' - + i^^j) + ^ Cmatter , (85) 

We see appearing the cosmological constant, the Einstein-Hilbert term with an extrinsic 
curvature term quite similar to the finite width corrections one obtains for cosmic strings 
in flat spacetime [67]! In other words, the only bulk quantity entering in the equations of 
motion is the topological defect fixing the overall mass scale and the extrinsic curvature of the 
surface giving a matter like component in the action. Therefore the Gauss-Bonnet term quite 
naturally gives an induced gravity term on the brane for codimension 2. This does not mean 
that there is a localised 0-modc graviton on the brane location but much like in DGP one 
can expect a quasi-localised 0-mode or in other words ordinary 4 dimensional gravity up to 
some cross-over scale as in the DGP model [31]. It would seem that the topological quantity 
/? would be giving us in this case a crossover scale for 4 dimensional gravity to 6 dimensional 
gravity. Clearly perturbation theory of well-defined warped backgrounds is needed in order 
to answer this question. 

Equation (79) has another surprising property. Indeed, an important simplification takes 
place if we suppose that 

[p/2] + 1 < n , (86) 
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ie that the codimension of the brane is larger than its intrinsic dimension. In that case, using 
Gauss's equation, we have 

<.*=(A";r)A«;t....„ (S7) 

and the matching conditions (79) arc simply the induced Lovelock equations on the branc with 
no extrinsic curvature terms! Therefore, the action for a distributional 3-brane embedded in 
D = 8, 10, ...2d dimensions is exactly the Einstein-Hilbert plus cosmological constant action 

^g>=3,n>l) ^ _ ^)Area(52n-l) J^^(an + an+lR'""") + Cmatter , 

with Planck scale set by 

Mil = (1 - /3) Area(52„-i) a„+i . 

For a 4 or 5-brane, we will have in addition the Gauss-Bonnet term etc. In other words, if 
the codimension verifies (86), all the extrinsic curvature corrections drop out and there is a 
complete Lovelock reduction from the bulk to the induced Lovelock terms on the brane. All 
the degrees of freedom originating from the bulk at zero thickness level are exactly integrated 
out giving the most general classical equations of motion for the brane. As a consequence we 
have in particular energy-conservation on the brane (see also [96]). 

4 Applications to braneworlds 

Braneworlds offer interesting and direct applications to Lovelock theory. This is due to 
the fact that Lovelock theory is the most general metric theory that can enjoy well defined 
junction conditions There are two reasons for this. Firstly, the field equations bccing of 
second order Dirac matter distributions still allow a continuous metric across the junction 
surface and therefore contraints are imposed on geometric quantities such as the extrinsic 
curvature and the induced curvature of the braneworld. They will lead inevitably to a self 
gravitating equation of motion for the surface in codimension 1 and in interesting contstraints 
for codimension 2. Furthermore, in higher codimension [97] Einstein theory does not pocces 
the necessary richness in order to admit distributional sources, rather the finite thickness of 
the defect plays a neccesary role or one has to admit the presence of non-removable curvature 
singularities at the brane location. Lovelock theory seemingly admits higher codimension 
distributions [45] and simple examples of codimension 4 braneworlds are emerging [47]. Here 
we will concentate in turn to codimension 1 and codimension 2 braneworlds. 

4.1 Codimension 1 braneworlds and their effective 4d gravity 

In this section^^ we will consider the application of Lovelock theory to a 4 dimensional 
braneworld embedded in a 5 dimensional curved backround spacetime. We will take a flat 

single brane in adS ie a Randall-Sundrum type braneworld and its perturbations and we will 
then give the relevant cosmological evolution equations. We will also include the induced 

thank Stephen Davis for unpubhshed collaboration throughout this section which in part can be found 
in his papers [68] 
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gravity term so as to consider the most general 5 dimensional configuration with a brane 
boundary. For simplicity we will also consider Z2 symmetry although the assymetric cases 
present a number of interesting features [27]. Consider therefore the action, 



S = — d^x^ {(R + aG-2A) 
2 Jm ^ ' 

' J d'^xV^ (K + 2a( J - 2G''^Kab) 



+ ^ ^ d^'xV^ [(5R - 2t) - ^ d^Xy/^C^at - j d^xV^Cmat 

where the first line represents the bulk Lovelock theory, the second the Gibbons-Hawking and 
Myers boundary terms and the third the induced gravity and matter contributions. Notice 
that we have included a dimensionless number in front of the Einstein-Hilbert term since 
we will take its zero limit in order to see the resulting 4 dimensional brane gravity of a pure 
Gauss-Bonnet bulk term. Given the normal vector rta, the 5 dimensional projector 

Kb = 9ab - naUb (89) 

is the surface induced metric and 

Kab = h'^aVcUb (90) 

is the extrinsic curvature (strictly speaking the first and second fundamental forms of E). 
The caret denotes tensors constructed out of hab- [X] as before denotes the jump in X across 
the brane. 

The field equations in the bulk and on the brane are, 

(Gab + aHab + Agab = M'^Tab (91) 

where Hab is is the second order Lovelock tensor (25). The bulk and brane energy-momentum 
tensors are respectively Tab = '^SCmat/Sh'^ - habJO-mat and Sab = '^SCrnat/Sh"''' - habJC-mat- 
According to (76) [38] we have 

C [Kab - Khab] + a [Sjab - Jhab + 2P„ed6i^'''] - PGab = -M-^Sab + Thab (92) 

For an adS bulk and a flat brane situated at ^; = we use a Poincare patch and the solution 
reads, 

ds'^ = a^{z)r]f,^dx''dx'' + dz^ (93) 
where the warp factor is a{z) = e"^'^' and we have, 

A = -6A;2(C-2aA;2) , T = 2k{'iC, - 2ak^) (94) 

for the bare cosmological constant and t brane tension respectively. The effective curvature 
scale k is given by 2ak'^ = C T a/C^ + 2aA/3 from solving the field equations in the bulk and 
on the boundary. The upper sign corresponds to the Einstein branch. 

To obtain the effective four-dimensional gravity induced on the brane, we consider a 
general linear perturbation theory parametrised by ^ab{x, z), around the background solution 
(93), 

ds^ = a^{z){r]i^^ + jf,^)dx^dx'' + 2-i^^dx^'dz + {I + '^zz)dz^ (95) 
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We will deal with gauge in a moment. We consider in particular the perturbed four dimen- 
sional Ricci tensor corresponding to the linear perturbation ^^n, 

= \ (25"5(.7,.)a - °47m- - d^d,^) (96) 

and we define TZ = r]^^TZ^y and Q^i, = TZ^u — (l/2)??^jy'7^. These are the typical geomet- 
rical quantities encountered in perturbation theory of induced gravity terms such as DGP 
braneworlds [31]. 

Given our mirror symmetry bulk equations (for z > 0) are then 

(C - 2afc2) {(Gfc^'^T.M - ^)e''' - 3/c(4fe7,, + d,-f)] = ^T,, (97) 
(C - 2ak^) {{d^d^-fzu - □47./.)e'''' " ^kd^lzz - 9.(5'^7/.v " 5^7) } = -^T,^ (98) 

and 

(C - 2aA;2) { {20^, + (r?^,n4 - d^d,)j,, - 2{d, - 2k){ri^,d'^v^ - 5(^t;,))) e^^^ 

^ {dl - 4kd,){j^, - rj^,^) + 3krif,,{d, - 4k)^,,} = -^e^^'T^u (99) 

The higher order Gauss-Bonnet contribution is recognised easily by the appearence of the 
coupling constant a while the usual Einstein terms are identified by ^. Note first of all how 
the metric perturbation in the bulk are quasi-identical to those in Einstein theory apart from 
the overall coupling of a with the warp factor k multiplying the differential operator. If there 
were no warp factor, for example in a Minkowski bulk, the higher order Gauss-Bonnet term 
would have had no contribution in linear perturbation theory. A general remark we can make 
is that since the coupling constant a is dimensionful a ~ [length]'^ needs inevitably a bulk 
curvature scale to couple to. This can be related for example to a bare cosmological constant 
in the bulk or an effective cosmological constant on the brane [29]. Also note that given 
our couplings to matter in action (88) we need to have Q — 2ak^ > 0. We otherwise quite 
clearly have a bulk ghost. If we switch off Einstein gravity in the bulk C = and choose 
a < we have Einstein like perturbation theory. When C, = 2ak'^ the linear perturbation 
operator switches off and seemingly we approach a strongly coupled limit. This is the case of 
Chern-Simons gravity which is very interesting by itself and we invite the interested reader 
to consult [58]. In turn the junction conditions at the brane z = are 

[(C - 2afe2) {3krif,^j,, - 2r]^,^d^Va + 2d^^v^) - d,{jf,^ - v^^ul)}] + 2{P + 4aA;)g^^ = -^S^^ 

(100) 

Here note that the higher order Lovelock term give two contributions. Firstly a Neumann 
type boundary term in the sense that it involves the first derivative of the perturbation metric 
and is exactly the same as the usual Einstein perturbation. Secondly, an induced gravity type 

term, which is accompanied by the induced gravity coupling (3 and which yields an induced 
Einstein term on the brane. Note that indices in the above equations are raised/lowered with 
<^ and that T\ = T,, + e'^'^^^'^r^^. 

^°When spacetime is not conformally invariant as for example for a black hole then the perturbation operator 
picks up extra terms related to the background Weyl curvature that appears in the field equations for Lovelock 
theory for fc > 2! 
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Before solving the above equation let us first deal with gauge freedom. Any infinitesimal 
bulk transformation, ^ + gives the transformation law ^^b ~^ lab + ^Cdib which 
leaves the bulk equations invariant. Here take the infitesimal shift, 

z z + e(x, z) 

x^" ^x^' + ^^'{x,z) + d^'C{x,z) (101) 
(where d^^'^ = 0) and obtain the linear isometries, 

h^^y ^ h^y - a(z)(^; + d^C') - d^e (102) 

hyy —> hyy — 2^ (l03) 

hp^v h^y - d^^jy - dyi^ - 2df,dy^ - 2a'{z)er]^y (104) 

which leave the bulk equations invariant. A bulk coordinate transformation however can also 
displace the brane by, F — > F + e(x, 2; = 0) where we call F the brane bending mode. This is 
avoided as long as we choose e{x, 2; = 0) = on the brane. 

It is easy now to solve the equations (97) and (98) when Tab = and /c 7^ by going to a 
fixed wall gauge where there is no brane bending F and the wall is maintained at z = 0. In 
this gauge the brane boundary conditions are also invariant under the bulk transformations 
[30]. We get, 

lzz = -^d,^ (105) 

l., = -^d,^ + B, (106) 

9"%, = (107) 

where 7^1/ = 7^1/ — (1/4)717^1/ is the trace-free part of 7^1/ and the trace reads 7 = r]^^i^y. We 
also have d^B^^ = and = but we can choose gauge for which = as we will do 

from now. 

The remaining bulk field equation is given by 

(C - 2ae) (dl - 4kd, + e'^^'D^^ = (108) 
for 2; > 0. The boundary conditions at z = give 

2(C - 2ak'')d,^^, + (/3 + 4aA:)n47/.. = " { "^/^^ " ^ (^m- - ^) ^} (109) 

and 

Ah 

(C + /3fe + 2aA;2)n47 = ^5 (110) 

It becomes clear here that, as we mentioned above, the equivalent of 'brane-bending' effects 

are included in the tracefull part of the metric 7 which is a genuine scalar mode since it 
couples to matter. Following Davis [38] (see also [29], [35]) we obtain the bulk solution (for 
> Q) which vanishes as z — 00, 

^^(p,.)~.^'W/.Jl^) (111) 
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for A; > 0, where K2 are special Bessel functions of the second kind. Therefore the boundary 
conditions imply 



K2ipe''\'\/k) 2k 1 



K^ip/k)' ' p^F{p/k)M3[^'"' sl"^^" ^p^'J-^l (112) 
where 

F{p/k) = 2(C - 2afe^)^^|^^ + {kP + 4ak') (113) 

and 

k 4 

^ " ~P^3(C + f5k + 2aP)M3'^ ^^^^^ 

Any zeros of F with respect to momentum will correspond to tachyon modes on the brane (see 
[29] for the 2 brane case) and this is essentially due to the fact that we have mixed boundary 
conditions. The condition for a perturbatively stable theory is to have no ghosts (brane or 
bulk) and no tachyons, therefore, C, — 2ak'^ > and pk + Aak'^ > 0. 

We now keep A; fixed. For large distances and very small momenta p/k 

F{p/k) = iC + kp + 2ak^) + (in ^ + 1e) + 0{p''/k^) (115) 

where ~ 0.577 is Euler's constant. For small distances and very large momenta p/k we 
have 

Fip/k) = {pk + 4ak^) + (C - 2afe2) I _ + 0{{p/k)-^) (116) 



Hence for large distances (small p) 



p p^ 



we obtain Einstein gravity with Afpj = M^(C + + 2ak'^)/k. This is true even if C = 
For the short distance behavior and large p, we define an effective scalar mode 

then we obtain 

Q,. - 2(r?,.n4 - d,d.)<i> = j^^r^^S,. + 0{p-') (119) 

□ 2(C-2aA;^) 1 

+ Pk + 2ak'^){P + Aak)NP ^ ' 

which is linearised Brans-Dicke gravity with Mpj = M^(/3 + 8a/c) and 

(--)-^%^ (-) 

oj a Brans-Dicke coupling. This is clearly an effect absent in conventinal GR and is present 
due to the Gauss-Bonnet term and the induced gravity term on the brane parametrised by 
a and /3 respectively. Solar system constraints dictate that a; > 4 10^ in order to agree with 
time delay experiments as that of the Cassini spacecraft, [75]. Therefore u should be pretty 
big and hence we need to have C, — 2ak'^ « 0, close to the Chern-Simons case. 



and 
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4.2 Codimension 1, brane cosmology 

Having looked at the perturbation theory of ra codimension 1 braneworld let us now study the 
exact solution for cosmology (see also Gregory's lecture notes [28]). Consider a 4-dimensional 
cosmological 3-brane of induced geometry, 

dsL = -dr^ + R\t) + X^d^i^ (122) 

We suppose, following the symmetries of our metric, that matter on the branc is modelled 
by a perfect fluid of energy density £ and pressure P. The brane is fixed at z = 0, and the 
energy-momentum tensor associated with the brane takes the form, 

5f)=diag(-^:(r),P(r),P(T),P(r)). 

We assume Z2 symmetry across the location of the brane at z = 0, and set = 1 for the 
time being. The bulk symmetries are 4 dimensional cosmological symmetries imposed from 
(122). In other words we have three dimensional isotropy and homogeneity which leads us 
to the bulk metric (26). Using (26) and solving the field equations in the bulk we found 
two hulk gauge degrees of freedom, U and V . These, leave the brane junction conditions 
invariant, or equivelantly, the distributional part of the field equations (25) invariant, if and 
only MU = V. The one remaining bulk-brane physical degree of freedom can be traced after 
coordinate transformation to the expansion factor or branc trajectory R = R{t) [71], [37] 
evolving in the black hole bulk (39). At the end of the day setting, 

where H = ^ the generalised Friedmann equation reads. 



f - 3M = \l ^"^""^^ ^ 2(4.4a + 2C + U) (123) 



where we have defined 



2 ^ 2ak 4a/x 



and we set 

f/o = J7(^ = 0) = C - ^ak^ (124) 

Thus to avoid bulk ghosts we must have [/q > 0, and so the lower branch is ruled out. The 
standard conservation equation on the brane remains valid and reads, 

— + 3iJ(P + f) = (125) 
ar 

Squaring (123) gives a third order polynomial for A, 

£ = 32a^3 + ^^'(32aC - 3/?') - ^(2C/2 -8(3- 8C')^ ' ' OiU + . (126) 

Setting £ = T + p and )U = we can read off the critical tension from (126). Critical 
means that the effective cosmological constant on the brane is zero, 

= 2 ^~^° (2C + Uof = 4fc2(3C - 2ak^f (127) 
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This condition has to be imposed if we want to analyse genuine geometric effect on the mod- 
ified Friedmann equation. In turn this gives to hnear order in p the reduced four-dimensional 
Planck mass 

' ' (128) 



^{2C - Uo) + PVC^ 
where we restore the fundamental mass scale M of the 5-dimensional theory. This agrees 
with (117) obtained in the previous section. 

Before proceeding to analysing specific cases it is useful to recast the Friedmann equation 
inputing Tc and the warp factor k. In order to avoid ghosts and tachyons we have (3 > —8ak/3 
so let us set (3 = Po — 8ak/3. Then (123) for k = reads, 

£ = SPoH^ - 8afc(l - Q)H^ + {l-M)^Q + ^Tc{l + ^) (129) 

where 



and 



i?4(§ - 4^)2 

Now we see that in accord with perturbation theory the Gauss-Bonnet term yields an ordinary 
Friedmann term just as in induced gravity. This is parametrised by /3o. We emphasize that 
this term is due to the higher order Lovelock correction and not due to the higher order 
Einstein-Hilbert term. Therefore a naive expectation that an Einstein-Hilbert term of higher 
dimension gives 4 dimensional gravity whereas a Gauss-Bonnet term differing phenomenology 
is not true. Quite the contrary, higher order Lovelock terms even in codimension 1 give 
naturally ordinary 4 d gravity at some scales. This can be also expected by the presence of 
the induced gravity term in the Myers boundary for Gauss-Bonnet theory . Starting from 
the above we can obtain the second FRW equation which will tell us about acceleration. 

As an example let us now consider the case of zero tension. For a general analysis see 
also [88]. Although this means that we fine tune the couplings, lak^ = 3^ in essence nothing 
special happens and the equations become easier to deal with. This means that in order not 
to have ghosts we take a < 0, upon which /?o ~ (for M=l). The Friedmann equation 
(for «; = 0) simplifies to, 

p = 3poH^ + 8|a|fc(l - Q)H^ - (1 - M)^^^Q (130) 

where now Q = \J ^ + ^^'x^ ^"^^ ^ y^l - -^^^^ with = 1 for a pure adS back- 
ground. For late time acceleration we need the second FRW equation. To obtain it, we 
differentiate (130) and successive use of (125) and (130) gives, 

8\a\k^ (H"^ , l-M\( , 2{1-M^)\ 16(1 - 7W)Q , , , 3 



(131) 



4 



P SlalA;^,, A\a\k f H"^ l-M 



thank Nemanja Kaloper for pointing this out 
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Setting a = gives ordinary Friedmann equations. This limit effectively kills all the 5 
dimensional part of the action for the special case we are treating. The second term in (130) 
can be interpreted as an a correction (of negative sign) and the third term a bulk black hole 
term (of positive sign). The terms on the LHS of (131) tell us whether there is possible 
acceleration or not at late time for ordinary matter {w = 0) or radiation w = 1/S. This fact 
is true modulo the sign of the parenthesis of the RHS which may also change sign as, we will 
see, inversing the accelerating equations of state but breaking the strong energy condition 
given (125). 

It is now indicative to set /U = and study late-time effects. Then (130) simplifies to, 



p = SPoH^ + l6\a\H^{k - ^ + H.'^) (132) 

and (131) to, 

p r p 4|a|/^2 ] ^^^^^ 



The first term on the RHS of (132) is the usual FRW term (with the right Planck mass) 
but the second term is of negative sign which means that we cannot take H arbitrarily large 
For late times cosmology approaches usual LFRW cosmology. Lastly let us set the bulk 
Einstein-Hilbert term to zero C = 0) take a < 0. Then Tc = —Aak^ and 



with 



M = ^|l+ ^ 



Then the effective LFRW equations reduces to, 

S = 3/3oH^ -8ak{l- Q)H^ -Aak^MQ (134) 

which gives at late times ordinary 4 dimensional cosmology despite the fact that there is no 
Einstein-Hilbert in the bulk action. This proves our claim in the introduction that the higher 
order terms of Lovelock theory enhance ordinary 4 dimensional gravity. 



4.3 Codimension 2 braneworlds 

Let us now look at a braneworld of codimension 2, in other words let us consider 6 dimen- 
sional bulk spacetime with 4 dimensional maximally symmetric subsections. The general bulk 
solutions, as we saw in section 2, are solitons of manifest axial symmetry with d/d6 as the 
angular Killing vector, 

ds^ = Vir)de^ + + r'^d^Ki (135) 

and d'^K/i, is the 4 dimensional line element of adS, flat or dS spacetime, n = —1, 0, 1 respec- 
tively. For suitable parameters in (57) the radial coordinate varies inbetween r_ < r < r-(-, 
where r = r± are the former horizon positions for (46) and will be the possible brane loca- 
tions, r = r± for (135). In particular r = r+ can be infinity itself in which case the effective 
volume element in the (r, 9) direction is infinite (for the analysis in the Einstein case see [72]). 
This is the codimension 2 version of warped compactification. We see therefore that the six 
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dimensional soliton (135) posseses the correct spacetime symmetries to describe the general 
maximally symmetric 4 dimensional braneworld of constant curvature. The Wick rotated 
version of the staticity theorem tells us in particular that axial symmetry comes for free and 
need not be imposed when resolving the system of bulk equations. In order to introduce 
codimension 2 branes carrying some 2 dimensional Dirac charge or tension we need to rein- 
troduce conical singularities at the relevant axis origins (59) which are at r = = r_ and 
r = Ta = r+ by allowing for the presence of deficit angles P±. Indeed from (60) we have the 
identification 0^0+ "^^^^^ which means that we have the relation, 



2^+ 2^- 



(136) 



which relates the topological parameters /3± with the geometrical quantities such as mass, 
charge of the soliton metric (135). In particular note that wc can always get rid of one 
of the conical singularities (and the resulting tensionful brane) and thus construct warped 
spacetimes with finite volume element and a single brane. Since the extrinsic curvatures for 
the branes are zero the brane junction conditions [44] are given by, 

27r(l - /3±) (s;: + 4aGl^^) = S; (137) 

with T^™"e _ Sfj^^6^'^\p) = S^i,^i^. We see the appearence of the induced Einstein tensor 
on the brane originating from the Gauss-Bonnet bulk term in the Lovelock action [45] . Note 
that the warp factor is given by the value of r^. and in particular for k = 1 say G^"*^ = 
—SH^j^i, = —-p-^^u where is the de sitter induced metric with curvature set to 1. The 
induced Newton's constant on the brane is given by 

A complete analysis of these solutions is given in [46] where self-tuning and self accelerating 
solutions are studied. 



5 The extended Kaluza-Klein reduction 

In the previous section we discussed warped compactifications of higher dimensional space- 
times. In this section we will discuss the case of Kaluza-Klcin cosmpactification. It is well 
known that the Kaluza-Klein (KK) reduction of Einstein theory gives us an Einstein-Maxwell- 
dilaton theory (EMD) for specific KK couplings and some periodic boundary conditions. But 
what is the resulting theory for the KK reduction of a higher dimensional Lovelock theory? 
Is the resulting reduced theory going to give second order field equations? Is the resulting 
theory unique as its higher dimensional counterpart? Let us take here for simplicity the 
Kaluza-Klein reduction of an Einstein-Gauss-Bonnet theory from D = 5 dimensions to d = 4 
dimensions. The full analysis in arbitrary D and for Lovelock theory was carried out in [48] 
(see also [14] for applications to cosmology). Starting with the relevant Lovelock theory in 5 
dimensions, 

S = — i- / (fx dy J-g(^) ( - 2A + aG^^A (139) 
IottG Jm ^ ^ ' 
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we consider the following Anzatz for the 5 dimensional metric, 

ds'^ = g'f^jdx'^dx'' 

= {9t.u + e-^'^Af,A^)dx''dx'' + 2Af,e-^'f'dx''dy + e'^^'^dy'^ 



(140) 



As usual we are making the basic assumption that there is a Killing vector dy in the 5th 
direction, in other words we do not consider here warped solutions we are rather interested in 
integrating out directly and obtaing the resulting 4 dimensional theory. We impose periodic 
boundary conditions on y. Wc expect that the resulting 4 dimensional theory will be an ex- 
tension of the Einstein Maxwell dilaton (EMD) theory in 4 dimensions with some exponential 
potential (in the presence of a 5 dimensional cosmological constant). Indeed integrating out 
the y direction we obtain [48], 

-^eff = / t^^^V^e"^^|i^-(V</')2-2Ae^'^-iF2 + aG(4) + 



+ 



IM4 

3a 

Te' 



2f;:f^f^fi^ 



int 



where note the presence of a non trivial interaction term which reads. 



'int 



= — I J d'^x^/—ge 



4F^,F-«i?t - F^R) 



(141) 



(142) 



The interesting result that one can prove is that the field equations obtained from variation of 
the fields are still of second order as for the higher dimensional Lovelock metric theory. In fact 
the higher order EMD theory in question is the most general second order theory that has up 
to second order partial derivatives. Any other numerical combination of the interaction terms 
for example would have given higher order derivatives. Hence we come to the interesting 
conclusion that higher dimensional Lovelock theories when dimensionally reduced via the KK 
formalism retain their nice properties dictated by Lovelock's theorem. In order to give the 
simple basic properties of the 4 dimensional theories we freeze out in turn the degrees of 
freedom in (141). 

Taking a constant dilaton the Gauss-Bonnet term drops out, since in 4 dimensions it is a 
topological invariant, and we are left with a modified Einstein-Maxwell theory, 

Seff = [ d^xV^{R-2A- + fl \{F^,Fn' - 2F;F,^F«F,H } - S^nt (143) 

Black hole solutions to this modified Einstein-Maxwell theory and for A = have been 
studied (partially numerically) [84] and they are corrected Reissner-Nordstrom solutions. For 
flat spacetime in particular, ie setting g^j^ = rj^i^, we get a non-linear version of Maxwell's 
theory which reads (see the nice paper by Kerner [85] whose notation we follow here), 



/ 



d^x- 



3a 



iF^^.F^"'f - 2F;F^^i^XF«^] } (144) 



with field equations {a = ■^^), 



0, 




(145) 
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We see that the usual Maxwell equations are corrected by non-linear terms with coupling 7. 
It is interesting to recast everything with respect to the electric and magnetic field, E, B. We 
then get, 



divE 
rot{B) 



3-7 ^ - ^ 
—B.gradiE 



dE_ 37 



B 



B) 

djE . B) 

di 



E A grad{E . B) 



div{B) = 
rot{E) = 



dB 

' dt 



(146) 



In this form we can make two obvious remarks. First of all whenever the electric and magnetic 
field are perpendicular to each other higher order terms drop out and hence usual EM solutions 
are unchanged. This holds in particular for electromagnetic wave solutions. However note 
that since we loose linearity one can no longer necessarilly superimpose electromagnetic wave 
solutions if they are not perpendicular to each other. Furthermore, we can define an induced 
charge density and current. 
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Jind — 2 



37 - - ^ ^ 

Pind = n-B . grad{E . B) 

B^^^^-EAgrad{E.B) 



(147) 



that simulate the higher order terms and verify the continuity equation, ^^f^ ~^div{jind) 



0. 



Let us now in turn freeze the vector field strength Fn,^ = 0. Then we obtain a second 



order scalar-tensor theory that has been studied quite a lot recently [38]. Numerical black 
hole solutions to such theories were discussed early on in [73] . It is important to note that in 
this case the Gauss-Bonnet 4-dimcnsional scalar is no longer redundant for the field equations 
and plays the role of the mediator in-between Einstein gravity and the scalar sector. This is 
true in whatever frame we choose to go. For example in the Einstein frame, in other words 
even when the scalar does not couple to linear order with gravity we have that it does so with 
the Gauss-Bonnet term. 



eff 



Mi 



{i?-(V</.)2 



(148) 



where d, 7 are some specific couplings depending on the dimensional reduction. Such models 

are higher order corrected dark energy-quintessence models that are used to evoke late-time 
acceleration of the universe. The essential point here is that although to leading order they do 
not affect solar system constaints once one includes the Gauss-Bonnet term this is no longer 
true and their solar system constraints can be rather stringent [86], [87]. Most stringent 
constraints arise from light time delay which is calculated from the Cassini spacecraft [75] . For 
cosmological constraints which are far weaker one can consult [87]. An interesting alternative 
has been put forward recently [83] in the case of a higher order generalisation of Brans-Dicke 
theory. It was there shown that the combined effect of the higher order corrections and the 
scalar sector can reduce or even eliminate the constraint on the Brans-Dicke parameter lo by 
imposing particular higher order coupling functions. 
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6 Concluding remarks and open problems 



In these short lectures we saw some of the basic properties of Lovelock theory. Some cosmo- 
logical applications were given for codimension 1 and 2 braneworlds. There are a large number 
of open problems in these theories and certain results /issues which we have not treated here. 
One of the aims of this closing section is to list and comment on some of these open problems. 

Let us start with exact solutions (check out the lecture notes of N Obers for black holes in 
Einstein gravity). For definiteness we studied only the static ones here, a Taub-NUT version of 
these has been found in [20] . Stationary metrics (akin to the Kerr black hole for 4 dimensional 
Einstein gravity) have not been found despite efforts. Another important solution which is 
missing is that of the black string (see [95] for a perturbative treatment). Although this 
solution is trivially found in GR in Lovelock theory this is not the case. The reason is that 
Ricci flat solutions in 4 dimensions (as is the 4 dimensional black hole in the 5 dimensional 
black string) are not vacuum solutions to Lovelock's equations. In fact the bulk Weyl tensor 
also contributes to the field equations (unlike in Einstein's equations!) for Lovelock order 
k >2 [88], [90] ie, once we switch on the Gauss-Bonnet invariant. Therefore a 4 dimensional 
Ricci flat solution which is non-conformally flat will not solve the Lovelock equations (see the 
nice analysis of [22]). The absence of such simple solutions is maybe an accident but maybe 
it is also questioning the relevance of the black string type of solutions which are already 
questionable [91]... At the same time if such an exact solution was found for Lovelock theory 
it would without doubt be genuinely different from its Einstein version. What could we then 
say of its stability? Another generic technical problem is that even in the absence of a bare 
cosmological constant in the action we always expect one in the Lovelock solutions due to the 
presence of multiple branches which have no flat Einstein theory limit! Thus the candidate 
metric for resolution must be written in an anzatz suitable for a cosmological constant solution 
even in the absence of a bare cosmological constant! We know from studies in Einstein theory 
that beyond a certain symmetry a cosmological constant spoils integrability [92], [93]. One 
way around this is to start by looking for solutions at special cases as the Chern-Simons case. 
Furthermore additional technical problems, due to the absence of Ricci flatness can arise if 
we simply translate Einstein results. For example, the diagonal Weyl anzatz in [94], [92] is a 
priori no longer true for k > 2 for in order to obtain it we use the fact that the background 
is Ricci flat^^! 

Another subject that deserves further attention is that of higher codimension braneworlds. 
We saw that in the context of Lovelock theory one could in principle define higher codimension 
braneworlds with Dirac distributions. Only recently [47] did we obtain the first example of 
a codimension 4 defect having only a removable Dirac singularity. In particular the case of 
codimension 2 and its cosmology has yet to be elucidated. In this review we explained how 
one could obtain the exact solutions describing the maximally symmetric branes [46] . To what 
extend does Lovelock theory permit us to recover Einstein gravity on the brane [44]? Our 
understanding from [45] is that Lovelock gravity gives us induced Einstein terms on the brane 
but not Einstein gravity on the brane at all scales. In other words similarily to DGP, only up 
to some crossover scale, do we expect gravity on a codimension 2 brane to be 4-dimcnsional. 
Localised 4 dimensional gravity will occur only when there exists a localised 4-dimensional 
zero mode graviton in the gravitational spectrum of perturbations. Therefore clearly, what 
is missing is a clear cut way of developping perturbation theory in Lovelock gravity. We 

^^I thank Robin Zegers for pointing this out 
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firmly believe that beyond the Gauss-Bonnet term one has to use differential form formalism 
which wc highlighted here. Only then will we know the gravitational spectrum of higher 
codimension braneworlds and their 4 dimensional phenomenology. We must point out that 
perturbations of Einstein-Gauss-Bonnet black holes has been carried out in [74] and causality 
issues have emerged in the context of the adS/GFT correspondance and EGB planar black 
hole perturbations [98]. 

We hope that with this manuscript we have communicated the fact that Lovelock theory 
is a technically challenging, interesting, well-motivated and exciting subject of research with 
numerous open problems that await resolution. 
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